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Abstract-- Kelvin effect (Skin effect) is used in surface 
hardening produced by induction heating of gears, 
camforms, camshafts and other work pieces of complex 
geometries. The induction heating equipment for surface 
hardening of metals and alloys (using LF or medium 
frequencies in the jargon of induction heating) is composed 
by a coil or coil assembly and a power semiconductor 
driving system up to 50kHz. The load seen by the driving 
system is equivalent to a transformer. The primary 
corresponds to the excitation coil or coil assembly, and the 
work piece corresponds to a short-circuited secondary. To 
asses the electrical load it is necessary to determine the 
variations in skin depth from place to place due to local 
curvature effects in the work piece, and its variations in 
space and time due to variations in conductivity and 
magnetic properties coupled with thermal effects. In these 
and others technical applications of Kelvin effect it is often 
necessary to be able to relate local skin depths with local 
curvatures of the surface of electrically conductive bodies. 
The purpose of this technical review is to show how to 
obtain a closed form analytical formula that relates the 
local skin depth with the local mean curvature and the 
well-known skin depth for a flat conductive body. Then, 
we suggest how to apply the above-mentioned formula to 
describe the electrical load of the induction heating 
equipment in the conditions used for surface hardening. In 
the choice or design of an induction heating system the 
parameters of the intended process (depth of Kelvin effect, 
temperatures to be reached and duration of the heating 
process, amongst others) put restrictions over the coils and 
the power driving system. To determine the best shape and 
size of induction coils or coil assemblies, the complex 
thermal and electromagnetic processes produced in the 
work pieces must be simulated using a computer code. The 
application of the formula given in this report, during the 
planning of these digital simulations is suggested. In the 
design of the electronic power system, driving the coils 
with the work pieces inside, the load impedance must be 
specified. The formula reviewed here could be applied to 
asses some characteristics of the load that may be of 
interest in the choice or design of the driving system. 


Index Terms-- Induction heating, Load modeling, Skin 
effect. 


(1)-INTRODUCTION 


Induction heating is a non-contact heating method. An 
electrically conductive work piece (a metal or a 
semiconductor) is located in an alternating magnetic 
field of a coil or coil assembly (Fig. 1). 


Fig.1 A simplified representation of the coil and the 
work piece. The total currents in the coil tube and in 
the workpiece are sketched. 


By electromagnetic induction the field produces eddy — 
currents in the material, which is heated because of Joule 
effect. 

In ferrous metals, hysteresis loss heating occurs, but this 
effect is in general fairly small compared with the 
heating effect of the eddy- currents. Hysteresis losses are 
considered elsewhere (See note 2 at the end of this 
report). 

The induced current density and its heating effects are 
non uniform. The current density is a maximum at the 
surface of the work piece nearest to the coil conductors. 
In the case of a simply connected body, the current 
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density is zero at its centre. 

As the local density of thermal power is proportional to 
the square of the local electric current density, the 
heating effect is also non - uniform. 

If the frequency of the alternating magnetic field is high 
enough, the induced electric currents are developed in a 
thin layer adjacent to the above-mentioned surface of the 
work piece. 

This is the skin or Kelvin effect, which is used in 
engineering manufacturing industry for the hardening of 
bearing surfaces, for welding, for brazing, and similar 
heat treatment processes [1]. 

The depth of penetration of eddy — currents can be 
characterized by the so called “skin depth”. It is a 
function of the frequency of the field, as well as the 
physical properties and geometry of the work piece [1]. 
For surface hardening, the material is heated during a 
few seconds. In the earlier stages of the process, the heat 
is confined in a layer with the skin depth. 

The heated area is then cooled very fast, so that a surface 
hardening effect is obtained while the bulk of the 
material remains ductile. 

The distribution of heating in the body can be controlled 
adjusting the frequency of the alternating magnetic field. 
Frequencies in the range 50 Hz to 1 MHz are used for 
purposes, through heating at 
frequencies to surface heating at high frequencies. 


several from low 
In this report, surface hardening using power sources up 
to 50 kHz will be considered. 

These are the so-called medium frequency heating 
equipments. They are based on high — power thyristors. 
In some of these equipments, the rectified 50 Hz supply 
is chopped by the thyristors and fed into a resonant 
circuit formed by the reactances of the loaded work coil 
(resistive and inductive) together with tuning capacitors. 
The load circuit seen by the driving system is equivalent 
to a short — circuited transformer. 

The primary corresponds to the excitation coil or coil 
assembly. The work piece can be considered as a short — 
circuited secondary. 

The impedance of this short circuited transformer 
depends on the geometry and physical properties of the 
work piece, the frequency of operation, and_ the 
geometry and number of turns of the coils. 

In the choice or design of an induction heating system to 
be applied in surface hardening processes, the 
parameters of the intended process (depth of Kelvin 
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effect, temperatures to be reached and duration of the 
heating process, amongst others) put restrictions over the 
coils and the driving system. 

In surface hardening and other technical applications of 
Kelvin effect, it is often necessary to be able to relate the 
variations of the skin depth from place to place, due to 
local curvature effects in the heated body, and its 
variations in space and time due to variations in 
conductivity and magnetic properties coupled with 
thermal effects. 

A complete description of the problems posed by the 
method can be found in reference [2]. 

The practical importance of finding a relation between 
the local skin depth and the local curvature of the 
boundary of the work piece was stressed in reference [3], 
where an interesting qualitative description is given. 

The purpose of this review is twofold. First, to show the 
mathematical deduction of closed form analytical 
formula that relates the local skin depth with the local 
mean curvature and with the well-known skin depth for a 
flat conductive body. Second, to suggest how to apply 
the above-mentioned formula to describe the load 
impedance that must be driven by the induction heating 
equipment in the usual conditions for surface hardening 
of gears and other work pieces of fairly complex 
geometries. 


(2)-AN ANALYTICAL FORMULA THAT RELATES SKIN 
DEPTH WITH LOCAL MEAN CURVATURE 


To study the Kelvin effect in low frequency fields, as in 
the case of induction heating, displacement current may 


be neglected [4]. In this case if S =E AH is the 
Poynting vector, from Maxwell equations we obtain': 

es ee es 

Ot 

Let us consider the work piece located inside the 

working coil. If an alternating magnetic field is produced 


() 


in the coil, the electromagnetic field in the air outside the 
work piece but inside the coil may be described by quasi 
- static equations. In the surface of the work piece the 


' In this technical report the vector product is represented 


by the symbol / (for example E \ H ) and the scalar 
product is represented by the symbol - (for example in 


the divergence V - SorinE-J ). 
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electromagnetic energy enters to the material along the 
normal to the considered point, due to the big difference 
in phase velocity between the air and the conductive 
body (refraction with absorption in the boundary) [5]. 

As consequence, the Poynting vector in the boundary of 
the work piece is perpendicular to this surface and points 
towards the interior. 


Let now us consider the lines of the field S inside the 
work piece, and let us suppose that there is, by 
continuity, at least locally, a congruence of surfaces 
orthogonal to these lines (one of these surfaces is the 
boundary of the heated body). See Fig. 2. 


Poynting field lines 


Fig.2. Vector lines corresponding to Poynting field 
and the associate congruence of surfaces near the 
boundary of the work piece. 


But S=S é, being és a unit vector and S the norm 


of the Poynting vector. One surface belonging to the 
above-mentioned congruence pass through the point in 


which 5S is located, and és is the normal vector to the 


surface at this point. According to a result of differential 
geometry and classical field theory applicable to a 
congruence of regular surfaces, the divergence of the 
normal vector is twice the mean curvature of the surface 
2 (6) (71. 

Then 


dS 


V-S =é,-VS+SV-é, =F, +259 (2) 
S 


Here —— is the directional derivative of S in the 


S 
direction of € gs» and s represents the arc length along the 


corresponding field line. From (1) and (2), substituting 
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J=oE and taking time averages it follows: 


OY 5g, ee ae 
ds Ot 


where o is the conductivity of the material. If the 


(3) 


eee OB 
magnetic hysteresis is negligible, H i may be put 
t 


equal to zero. Then it follows that: 


d(n5) oe? 
= 260 a (4) 


The spatial scales that characterize the variation of the 
fields in a direction tangent to the boundary of the work 
piece are of the same order of magnitude of a 
characteristic dimension of the heated body. The spatial 
scale of variation of the fields along the direction of the 
Poynting vectors is of the same order of magnitude of 
the skin depth 6 for a flat boundary [4]: 


(5) 
OU 

Here @ is the angular frequency of the alternating 

magnetic field produced by the equipment and u is the 

magnetic permeability of the material. 

In the case of surface hardening, the skin depth 6 is 
several orders of magnitude less than a characteristic 
dimension of the heated body. Therefore, in a first 
approximation, the electric and magnetic fields can be 
described by the formulae corresponding to a flat 
boundary [4] [8]. The following equations are verified: 


E? =(E)? and S =——(EVH). 


V2 


Here ( ) represents the RMS value of the magnitude 


1 
that appears between these brackets, and the factor a 


appears because the electric field has a phase lag of 1/4 
relative to the magnetic field (neglecting magnetic 
hysteresis). 

Between the electric and magnetic fields there is the 
relation ([4] [8]): 

(E)=¢(H) 6) 

The parameter ¢ is given as the following function of 
and the 


the angular frequency, the permeability, 


conductivity: 
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@ 
ic = aselas (7) 
oO 

From (4) to (7) it follows that: 

d(InS y) 
d(inS) =—20 -— (8) 

ds 
If the boundary is flat, (9 = O and from (8): 
a8, 

S(s)=S(0)e ® (9) 


Equation (8) gives the vanishing of the average 
magnitude of the Poynting vector going towards the 
interior of the body from its surface. The magnitude of 
the average Poynting vector in a point of the surface of 


the work piece is S (0), and S’ (s) is the corresponding 
magnitude at a point of the same field line at a distance s 
form the surface. This result suggests the introduction of 
an equivalent skin depth 6, that corresponds to the case 
of a curved surface: 


ra 
1+d¢ 


+ == or 6, = (10) 


5. 
This formula gives the relation between the local skin 
depth, the local mean curvature of the boundary and the 
skin depth corresponding to a flat boundary between the 
work piece and the surrounding air. 


(3)- THE LOAD IMPEDANCE SEEN BY THE INDUCTION 
HEATING POWER SOURCE 


Now, let us consider the work piece as a short-circuited 
secondary of a transformer equivalent to the loaded coil. 
In the conditions prevailing during surface hardening the 
skin depth is than the geometric 


dimensions of the heated body. In this case we can 


much smaller 
consider a bar of non-circular cross section. The 
boundary of each cross section is given by the same 
curve IT’ that may be described in the x-z plane of an 
orthogonal system of coordinates x,y,z. The y axis is 
parallel to the axis of the bar. Let h be the length of the 
bar. It is possible to introduce a surface density of 
electric current integrating the volumetric density in 
depth following the lines of the Poynting vector field. If 


Aa 


de 


s 18 the phasor that represents the surface density, and 


A 


[,,, is the global current that is induced in the work 
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piece, then: 
(1) 
The electric currents flow parallel to the x-z plane, in 
each cross section of the bar. 
If E 51S the phasor that represents the electric field 


tangent to the boundary of the heated body (also 
contained in the x-z plane), the induced electromotive 
force is given by: 


3=fE,dl (12) 


In the case of a half infinite medium with a flat 
boundary, it is possible to show that [9] [10]: 


é,-(E2), 
(Eb 


If the Kelvin effect is strong enough, and the curvature 


(13) 


2 is less than 1/5, we can suppose that (13) may be 


applied to the bar substituting 5 by the local skin depth 
de given by (10). Then multiplying and dividing by h and 
taking into account (11) it follows: 


ey eee oe ee 
a {eisah, 


So we obtain the following formula for the impedance of 


(14) 


the work piece: 


7 2, -[ S262) 


If I() is the length of the curve I’, from (10) and (15) 


we obtain 


(15) 


i, A) 
Lyp = — (1+ 7 (16) 
"oh, ( ) 
Here, by definition 
5 ae (17) 
1+é6¢0 


This 6, is an average skin depth that corresponds to an 


average mean curvature: 


— bedl (18) 
ir)! 


To estimate the electrical load seen by the induction 
heating equipment, it is necessary to calculate the 


impedance of the loaded coil. If a voltage V,. is applied 


Cc 


to the coil terminals and if /. is the electric current in 
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the coil, R, is the coil’s resistance, L, ,is the internal 


coil’s inductance, and @ is the magnetic flux through the 
coil’s interior ,and N is the number of turns of the coil, 
then we suppose that: 


I 
V. =Rpat 2 < Fe 
" dt dt 


(19) 


If®, 


interior surface of the coil and the surface of the work 


is the magnetic flux in the air space between the 


piece , and @® ,,,, is the magnetic flux that goes through 


the work piece, then: 


=0,+0,, (20) 


If A, is the area of the coil’s cross section and A... is the 


area of the work piece cross section: 


®, =A,(l-v)u,H 


a 


(21) 


oO 


wp 


Herev = is the so-called filling factor of the work 


[33 


piece relative to the coil interior. 
eee NI 
The mean field is H, = a being / the length of the 


system. 

Now, let us work with phasors. 

Then the complex electromotive force that is induced in 
the work piece is given by 


d®,, 
ar 


A 


(22) 


The electromotive force reflected into the coil is: 
GOs as 
= N d = 7 Ny 
t 


But between the current [ 


Cc 


in the coil and the 


current/,,, there is the transformer current relation: 


'P 


i,, =-NI -- Consequently, the fem reflected over the 


wp 


coil is: 
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d® , Aer 
-N—* =-Z,,N°I, (23) 
dt 
From (21) it follows that 
d® N° > 
N —- = jaA_ \(l-—v) —I 24 
jee) (24) 


From (19), (20) , (23) and (24) it follows: 
wp” c 


2 
V.=R I, + joLk,1.+ jo ~i, +N°Z 1 


Then, 
impedance, 


using equation (16) for 
the 
impedance Z, = R, + jX , the equations: 


) 


the work piece 


we derive for complete load 


R, =R.+ N°’ £ (25) 

oe o ho, 
Ne. Soe UE) 

X, =OL.; + J0A,U-v +N = (26) 
L ivJd ( ) Fs = he. 

As t+ ea, re : fedli 
Ss —_-=-2 an = ST 1S an average 

ro} }, ur)? 


mean curvature, both the active resistance R, and the 
reactance X, of the load are growing functions of the 


average curvature. 
Xx; 
The e inductance L, = —~ has three terms. 
@O 


The first is related with the self inductance of the coil, 
and it suffers the influence of the Kelvin effect. 

The second term is relatively independent of@ . 

The third term is the sum of two inductances. 

An inductance arising fromdé, that is proportional 


1 fs os 
to ——. Another inductance arising from the effect of 0 , 


Jo 


. 1 
that is proportional to— . 
(a) 


Both inductances decrease whenq@ increases, but the 
component related with the curvature decreases faster 
than the component related with the skin depth 6 of the 
flat body. 


The coil resistance R, also varies with the frequency of 
the alternating current due to the skin effect in the coil 


tubes. 
The effective resistance, as shown in (25), is the sum of 
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this resistance R., and a resistance numerically identical 


KC) 
o ho, 


This stems from the variation of the phase of the electric 


to the third term of the effective reactance N” 


currents with the variation in depth from the surface of 
the work piece, relative to the phase of the electric field 
at the surface. 


(4)-DISCUSSION AND CONCLUSIONS 


To apply the results obtained above for the skin depth 
and for the impedance of a work piece, we will consider 
first the simplest and more common type of gear: spur 
gears [11]. Spur gears are used to transmit rotary 
motions between parallel shafts, with gears teeth that are 
located on the circumference of a circular plate. At the 
teeth level the boundary is convex. But if the boundary is 
convex, the mean curvature is negative. Then, the local 
skin depth, according to formula (10) will be greater 
than the skin depth of a flat boundary. If the boundary is 
concave, the mean curvature will be positive. Then the 
local skin depth will be smaller than the skin depth of a 
flat boundary. This is shown in Fig. 3. In this figure the 
arrows that point towards the material of the work piece 
represent the field of Poynting vectors. The arrows, 
tangent to the boundary, represent the electric field. The 
magnetic field is perpendicular to the plane of the 
picture. 


Fig.3. Variations of the local skin depths with the 
local curvatures of the boundary for a work piece. 
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These variations in local skin depth are according to our 
intuition .If the adjacent boundary is convex, the field of 
Poynting vectors is convergent towards the interior and 
thus it should be less attenuated than the field adjacent to 
a flat boundary. If the boundary is concave, the field of 
power density flux will diverge towards the interior, and 
it should be attenuated more than the field adjacent to a 
flat boundary [3]. 

Now, let us suppose that this toothed wheel is located 
inside a coil to be heated by induction. The induced 
electric current will circulate in a layer adjacent to the 
surface of the teeth and the notches between the teeth, 
vanishing towards the interior of the wheel. As the local 
skin depth is higher at the teeth level than at the notches 
level, the real part of the local value of the impedance 
(the so-called active resistance) will be lower at the teeth 
level than at the notches level. However, the same 
electric current circulates everywhere (so that the electric 
current density at the teeth level will be smaller than the 
electric current density at the notches level). As a 
consequence, the material adjacent to the notches will be 
heated more than the material of the teeth. This fact is 
important from the standpoint of induction hardening. 
Moreover, from the results derived in this paper it is 
possible to calculate the distribution of heat power in a 
layer adjacent to the surface of the work piece. This 
distribution may be used as an input in the equation of 
heat transfer, and solving this equation the 
corresponding temperature field may be calculated. 

If the frequency of operation is high enough, the skin 
depth for a flat boundary 6 will be so small that dg 


will be negligible in comparison with 1. In that case, the 
local skin depth 6, will be equal to 6 and the heat 


power produced by the induced electric current will be 
the same in all locations of the surface of the gear. 

These consequences of the analytical formulae are 
confirmed both by experimental results and by digital 
simulations of the process of surface hardening of gears 
by induction heating. [2] [3]. 

The assumptions made to derive (10) suggest that the 
formula for the local skin depth can be applied only if 


the product 6 \o| is small enough, perhaps less than 0.3. 


This guess may be substantiated in a comparison 
between the approximate analytical formulae obtained in 
this paper for bars of any cross sections and the exact 
results obtained for the Kelvin effect in infinite circular 
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cylinders, in two different situations. 

In one of them, an infinite cylindrical conductor 
carries an alternating current parallel to its axis. The 
currents produce an field 
perpendicular to the axis of the cylinder [3] [8]. 

In the other situation, an infinite cylindrical conductor 


alternating magnetic 


is in an alternating external magnetic field, uniform and 
parallel to its axis. The induced currents are contained in 
the circular cross sections of the cylinder. 


This case, studied by Foérster and others [12], is a 
simplified mathematical model that may be applied in 
the framework of certain electromagnetic methods for 
nondestructive testing of metals. However, it has a closer 
connection with the induction heating case. 

A detailed analysis of the relation between the 
and the 
of the exact solutions for fields and 
currents in these cases will be done elsewhere. 


approximate formulae obtained here 


consequences 


Here we only note that for a convex cylinder the mean 
curvature is constant, is negative and in absolute value it 


is equal to half the reciprocal of the radius7, of the 


cylinder, so that from (10) it follows: 


6, = : 


1-6{ 1) 
2F, 


Exact solutions for the fields in 


(27) 


infinite cylinders, 
both the classical one for the Kelvin effect in an 
alternating current conductor and the solution obtained 
by Forster for the conductor inside an infinite coil, are 
constructed using Bessel functions [8] [9] [12]. 

From an analysis of the asymptotic expansions of this 


Nh r 
Bessel functions, if a and 5 are big relative to one, it 


follows that the amplitudes of the electric and magnetic 
fields inside the cylinder vanish toward the axis like 


(28) 


It is possible to make the approximation: 
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sl si ad 
to 2 1‘ =@ 2 'D ~ 


=") 
21% 


e (29) 


This relation is valid when the quotient between the 
depth from the surfacer, —r and the radiusr, is small 
enough relative to one. 


From Equation (20) and (21) we obtain that the fields 
vanish toward the axis of the cylinder like: 


=") _(»-r) em O (m-r) 
é 21% Z cS é fo) 21% = Op (30) 
Thus, we derive Equation. (24). 

Equation (27) may be used to estimate the way in which 
the amplitudes of the local fields and currents vanish 
with depth from the boundary of the work piece. 

In Fig. 4 several curves show the function: 


ed {5} 


=e (31) 
The absolute value of K is equal to 6/r0 .The different 
curves shown correspond to K=-1, -0.5, -0.2, 0, 0.2, 0.5, 
1. 


; 
: 
i 
: 
; 
5 
i 


1 15 2 25 3 


Dist to surface relative to skindepth (flat body) 


Fig.4. Electric current density at a certain depth from 
the surface, relative to its value on the surface as a 
function of the quotient between depth and the skin 
depth of a flat body. 


Positive values of K correspond to concave cylindrical 
surfaces and negative values correspond to convex 
cylindrical surfaces. An example may be seen in Fig. 5. 
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Fig.5 A coaxial cable’s cross section. The cable is 
formed by a solid internal cylinder and a hollow 
external cylinder. Mean curvature of the convex 
internal cylinder is negative and mean curvature of 
the concave external cylinder is positive. 


Fig. 6 shows the electric field inside the convex internal 
cylinder and the concave external cylinder for the 
coaxial conductor shown in Fig. 5. 


Fig. 6 The quotient between the electric field at a 
given depth and the electric field on the surface, as a 
function of the quotient between the depth from the 
surface and the skin depth for a flat body. 


In this case the sign convention is opposite to that 
used in Fig.4.Positive values of K correspond now to 
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concave external cylinders and negative values 
correspond now to convex internal cylinders. 

This figure represents the exact solutions given in terms 
of Bessel functions [3]. The exact results shown in Fig. 6 
may be compared with the approximate results given in 
Fig. 5, obtained using Equation (19) for the local skin 
depth. 

Formula (16) gives an analytical estimate of the work 
piece impedance. If the line integral of the local mean 
curvature, taken over the curvel’, is negative (as is the 
case of a predominantly convex body) both the resistive 
and the inductive parts of the impedance will be less 
than the corresponding impedance components for a flat 
body. If the line integral is positive (as is the case of a 
predominantly concave body) both the resistive and 
inductive parts will be greater than the corresponding 
impedance components of a flat body. If the frequency is 
high enough, these differences vanish. 

In order to calculate the impedance seen by the 
electronic driving system, the work piece impedance was 
reflected from the secondary to the primary circuit (the 
coil) of the equivalent transformer, using a well-known 
procedure. Equations (25) and (26) were thus derived. 
They correspond to the simplest mathematical model of 
the electric load. End effects and leakage of fluxes were 
neglected in this analytical approach, as well as the 
variations in conductivity and magnetic permeability 
produced by local temperature increase and magnetic 
hysteresis. 

Hysteresis may be described using a complex magnetic 
permeability in the formula for the skin depth [9]. 

One consequence is that the active resistance is now 
higher than the inductive reactance in the impedance 
given by (16). Magnetic saturation may be considered 
using a suitably smaller value for the magnetic 
permeability in the formulae for the skin depth and for 
the work piece impedance. 

Anyhow, to determine the best shape and size of 
induction coils or coil assemblies, in order to heat by 
induction, the work pieces of complex geometry as 
found in practice, the analytical approach is clearly 
insufficient. The complex thermal and electromagnetic 
and the 
interaction between the work pieces and the coils must 


processes produced in the work pieces 


be simulated using computer codes. The use of these 
codes is not straight forward, and wrong results can be 
obtained from the numerical calculations. 
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Nevertheless, the analytical approach affords general 
results that may be used as guidelines about what to 
expect and about what to search in relation with the 
results of the digital simulations. 


Also, the analytical approach gives us useful lower 
and upper bounds for the values of parameters such as 
the local skin depths or the impedances that may be used 
to select induction heat equipment and to establish some 
main characteristics of the heating process. All this 
seems to deserve further study. * 
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